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ABSTRACT 

The  Lindstrom-Madden  method  of  computing  lower  confidence  limits  for 
series  systems  with  unlike  components  is  extended  to  series  systems  with 
repeated  components  utilizing  the  results  of  Harris  and  Soms  (1983)^.  An  exact 
solution  is  given  for  no  failures  and  key  test  results,  together  with  an 
approximation  for  the  general  case.  Numerical  examples  are  also  provided. 
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significance  and  explanation 


series  systems  with  repeated  components  arise  often  in  engineering  and 
physics.  It  is  therefore  Important  to  utilize  data  obtained  on  individual 
components  in  an  efficient  manner  when  assessing  the  reliability  of  the 
■  .rrhined  system.  This  paper  gives  one  method  for  doing  so. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  Ires  with  MFC ,  and  not  with  the  author  of  this  report. 


THE  LINDSTROM-MADDEN  METHOD  FOR  SERIES  SYSTEMS  WITH  REPEATED  COMPONENTS 


Andrew  p.  Some* 


1.  INTRODUCTION  AND  SUMMARY 

A  problem  of  substantial  importance  to  practitioners  in  reliability  is  the 
statistical  estimation  of  the  reliability  of  a  series  system  of  stochastically  independent 
components  when  some  components  are  repeated,  using  experimental  data  collected  on  the 
individual  components.  In  the  situations  discussed  in  this  paper,  the  component  data 
consist  of  a  sequence  of  Bernoulli  trials.  Thus,  for  component  i,  i  ■  1,2, ...,lc,  the 
data  is  the  pair  where  n^^  is  the  number  of  trials  and  Y^  is  the  number  of 

observations  for  which  the  component  functions.  Yf ,Y2, • « • ,Y^  are  assumed  to  be  mutually 
independent  random  variables.  We  assume  that  there  are  Y^  components  of  type  i, 

1  <  i  <  k.  Then  the  parameter  of  interest  is  h(p1 ,p2, . . . ,p^)  “  h(p),  the  reliability  of 
the  system,  where 


h(p) 


More  specifically,  it  is  desired  to  obtain  a  Buehler  (1957)  optimal  lower  1  -  a 
confidence  limit  on  h(p). 

The  case  of  Y 1  =  Y2  =  •• •  =  Yk  “  1  has  been  treated  in  Sudakov  (1974),  Winterbottom 
(1974),  and  Harris  and  Soms  (1983). 

In  Section  2  we  summarize  the  general  theory  of  Harris  and  Soms  (1983)  applicable 
here.  In  Section  3  the  exact  solutions  to  no  failures  and  key  test  results  are  given. 
Lindstrom-Madden  type  approximations  are  given  in  Section  4.  Section  5  contains  numerical 
examples. 
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Bi'FHLEK'S  METHOD  FOR  OPTIMAL  CONFIDENCE  LIMITS 
We  now  specialize  the  general  results  of  Harris  and  Soms  (1983)  on  optimal  confidence 
limits  for  system  reliability  to  a  series  system  with  independent  and  repented 
components.  As  in  Section  1,  let 
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0  <  p,  <  1,  Xj  =  ni  ■  Yi(  xi  =  nj^  -  yi(  1  <  i  <  k,  S  =  {x^  *  0,1,...,^,  1  <  i  <  k}  and 
let  g(x)  =  ( x ! ,x2 , . . . ,x^)  be  an  ordering  function,  i.e.,  for  real  x^,  0  <  x^  <  n^, 
g(x)  is  non-decreasing  in  each  component.  It  is  often  convenient  to  normalize  g(x)  by 
letting  g ( 0 )  =  1  and  g(n)  =  0.  With  such  a  normalization,  g(x)  is  often  selected  to 
be  a  point  estimator  of  h(p).  Also  let  R  =  {r^,r  , .».,r  ,  s  >  2}  be  the  range  set  of 
g(x).  with  no  loss  of  generality  we  order  R  so  that  rj  >  T2  >  ...  >  rs  and  let 
A^  =  (x|g(x)  =  r^,  x  e  s,  i  =  1,2,...,s}.  The  sets  A^  constitute  a  partition  of  S 
induced  by  g(x).  we  assume  throughout  that  the  data  is  distributed  by 


n .  n .  -x .  x . 


f  ( x 


;p)  =  P~(X  =  X)  “  TT  1qi1 

y  i=1  i 
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(2.1  ) 


where  qi  =  1  -  pi(  i  =  1,2,...,k.  With  no  loss  of  generality,  we  assume 
n1  <  n2  <  ...  <  nk. 

From  these  definitions,  it  follows  that 


j 

P'jx  c  M  A.}  =  P~{g(X)  >  r.}  . 
P  iTl  1  p  3 


(2.2) 


From  (2.1)  and  (2.2),  we  have 
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P'tg(X)  >  r  }  *  l  l  •••  I  f(i;p)  ,  (2.3 

ir°  V0  \=° 

where  i  =  (i  i , i-2’  * "  * '  *k^  an<*  u2  =  u2^1^,',*,uk  =  uk  ^  1 '  ^2 '  •  •  • '  ^k- 1  ^  are  integers 
determined  by  rj.  Equivalently, 

[t1]  ft2]  [V 

P2(g(X)  >  r  }  =  l  l  •••  \  f(i,p)  ,  (2.4 

F  J  i ,  -0  i  =0  i,  =0 
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where  t2  =*  t2(  i tk  «  tj<(  i  1  ,i2  ,  •  •  • ,  ik_  i ) ,  with  t.|  =  sup{tl0  <  t  <  and 

g(  t  ,0 ,0 , . . .  ,0 )  >  r^}  and  t ^(i1 , i2 , . . . ,i 1 )  -  sup{t|0  <  t  <  n^  and 
g(i1,i2,...,ii_1,t,0,...,0)  >  r  },  l  =  2,3,.  ..,k. 

We  now  introduce  the  notion  of  Buehler  optimal  confidence  limits.  Let  g(x)  =  r 


Then  define 


a^^j  -  inf {h(p) | P~{i | g(i )  >  g(x)}  >  a} 


Equivalently,  by  (2.2),  we  can  also  write 


aqU)  *  infih<Pllp;ix  e  LJ  A.  }  >  a  1  . 

y  p  i=1 


Then  we  have,  from  Harris  and  Soms  (19S3), 


Theorem  2.1.  a  ,■*.  is  a  1  -  a  lower  confidence  limit  for  h(p).  if  b  ~  is  an' 
-  g(x)  g(X) 

other  1  -  a  lower  confidence  limit  for  h(p)  with  b  >  b  >  ...  >  b  ,  then 

1  r2  3 

b  <  a  for  all  x  e  s. 

g(x)  g(x) 

Two  possible  choices  of  g(x)  are 
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g(x)  =  tt  ft  c  \  *  U--J  . 
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Both  reduce  to  the  generally  used  g(x)  for  series  systems  with  independent  components 
when  Y1  =  Y2  =  ...  =  Yk=  1,  i.e.. 


g(x)  =  *1  [  tnA  -  XjJ/nj  . 

i=1 


Since  (2.7)  is  the  maximum  likelihood  estimator  of  h(p)  we  will  use  it  here  and  from  now 
on  it  will  be  understood  that  g(x)  is  given  by  (2.7).  With  this  choice  of  g(x),  we 
assume  from  now  on  that  0  <  x^  <  n^,  i  =  1,2,...,k,  since  ag(x)  =  0  if  some 
xi  =  ni*  With  this  assumption,  the  t±  in  (2.4)  are  given  by 


k  Y.  /  k  Y.  1/Y. 
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(2.10) 


t  =  2, . . . ,k,  with 


k  Y. 

FT  n.1  =  1. 

i=k+1 


For  the  purpose  of  simplifying  the  calculation  of  ag(J)  in  special  cases  it  is 
necessary  to  state  additional  results  from  Harris  and  Soms  (1983). 


Theorem  2.2.  Let  g(x)  ■  rj  and  let 


f*(x;a)  =  sup  P-{g(X)  >  r.},  0  <  a  <  1 

h(p)*a  P  3 


(2.11 ) 
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inf  f  (xi a)  *  0,  sup  f  (xja)  "  1 
0<a<1  0<a< 1 


and  f  (x;a)  is  strictly  increasing  in  a. 

*  •*. 

Theorem  2.3.  f  (xja)  =  a  has  exactly  one  solution  a  in  a  and  a  =  a 
-  1  a  a  g(x) 


3.  EXACT  SOLUTIONS  FOR  ZERO  FAILURES  AND  KEY  TEST  RESULTS 

We  first  assume  that  x  =  (0,0,.. .,0)  »  0  and  use  Theorem  2.3  to  obtain  a  ~ 

g(  o ) 

Theorem  3.1.  If  x  “  0,  then 
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f  (0;a)  =  sup  II  Pi  “ 


K  Yi  1-1 
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where  n,/Y.  =  min  n./y.  and 
3  3  1<i<k  1  1 


a  "•  ■  a 

g(0) 


i-1  i-1  i=1 


since  n^Y^  “  "j^i  *  0  *s  ecIuiva3-ent  to  n^/Y^  *  nj/Yj'  which  is  true,  and  therefore 

k  ( n^Y  ,-n  ,Y^)/Y .  1/Y  . 

1  T  p^  1  3  31  <1.  (3.1)  follows  by  noting  that  the  choice  p^  =  a  3,  p^^  =  1, 


i  *  j,  gives  |  |  p^1  =  a  .  Then,  using  Theorem  2.3,  we  obtain  (3.2),  which  reduces 

i=1 

to  the  known  series  result  if  Y1  =  Yj  *  • • •  “  Yk  “  1 • 


*-•  ■>  .v ^ > :>• 


where  g(x)  *  r ^ ,  by  the  monotone  likelihood  ratio  property  of  the  beta  distribution 
(Ix(a,b)  has  a  monotone  likelihood  ratio  in  -a  for  fixed  b,  which  implies  that 
Ix(a,b)  is  a  decreasing  function  of  a).  A  similar  argument  applies  to  the  other 
indexes.  Thus,  if  (3.4)  is  true  for  n,  =  n2  ~  •••  =  nk '  (3.5)  it  follows  for 

n,  <  n2  <  ...  <  nk. 

So,  assuming  n  *  (n, ,n, , . . . ,n,  ) ,  we  seek  to  maximize 
k  n1  k  k 

Ppi  l  1  * >  1  <"i  -  xi>  =  l  y*}  -  (3-6) 

p  i=1  j=1  J  i= 1  i=1 


where  Y^j  are  independent  Bernoulli  random  variables  with  parameter  p^  and 
k  Y .  k  y .  k  1/Y  . 

TTp  1  =  a.  If  TT  P, 1  =  a,  then  |  p,  ranges  from  a  3  to 


i*1 

1/Y 
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a  ,  Y^  =  min  Y^  This  is  seen  as  follows: 


1<i<k 


k  _k  y  1/Y,  k  1-yVy, 
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=  a  3  TT  Pi  3  1  3  >  a  3 

i=1 


i*  3 
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and  the  choices  p,  =  a  ,  p,  =  •  • .  =  P],  =  1 ,  and  p.i  =  a  »  Pi  =  1,  i  *  j»  attain 

k 

these  values.  From  the  results  of  Pledger  and  Proschan  (1971),  for  each  b  =  ]  f  pif 

i=1 
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a  <  b  <  a  ,  (3.6)  is  maximized  by  p1  =  b,  pi  =  1 ,  2  <  i  *  k.  Further,  the 

1/Yj  1/Y,  1/Y, 

maximum  over  b,  a  <  b  <  a  ,  of  the  maxima  for  each  b  is  given  by  p^  =  a  , 


pi=1,2<i<k,  by  the  monotone  likelihood  ratio  property  of  the  binomial 


1/Y, 


distribution,  and  p^  =  a 
the  proof. 

=  Yl 


p^  =  1,  2  <  i  <  k,  satisfies  |  |  p.  =  a.  This  completes 

i=1 


If  Y,  =  Y2  = 


,k  -  1,  some  guidelines  for  the  verification  of  (3.3)  are  given 
in  Harris  and  Soms  (1983).  In  the  present  case  (3.3)  must  be  verified  by  trial  and  error 

n 

i=1  ‘  i=2 


by  showing  that  min 
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IT  (nf  -  xf)  =  (n,  -  xf)  TT  nf 


and  that 
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max  |  |  (nA  -  x.)  <  (n1  -  x1 )  1  f  n.  . 
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1  Xi=X1+1 
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Example  3.1.  Let  k=  3,  n  =  (5,5,5),  Y  ”  (3,3,2),  a=  .10  and  x  =  (1,0,0).  Then 

3  rA  3  r, 

min  ^  |”  (n^  -  x.)  =  200000  and  max  ~\  f  (n^  -  x.)  «  140625  and 

■ .  .  *  ,  ■  . .  a 


3  i=1 

)  x.  =  1 

■  L.  a 


3  i-1 

i  *i=2 
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i=1  i=1 

x  is  a  key  test  result  and  (3.3)  is  satisfied  and  hence 


agU)  “  -41613  =  -0720  # 


where  .10  =  I.4161<4'2)-  Further,  it  can  also  be  verified  that  x  =  (2,0,0)  is  a  key 


test  result  for  which  (3.3)  is  satisfied,  but  that  for  x  =  (3,0,0),  (3.3)  is  violated. 


Note  that  Theorem  3.2  asserts  that  a  ~  =  b  for  0  <  a  <  1.  It  is  thus  possible 

g(x)  a 

that  (3.3)  is  not  true  but  the  conclusion  still  holds  for  o  of  practical  importance. 

This  is  taken  up  in  Section  4. 


4.  THE  LINDSTROM-MADDEN  METHOD  FOR  SERIES  SYSTEMS  WITH 
REPEATED  COMPONENTS 


When  Y1=Y2=***=Yr=1,  the  Lindstrom-Madden  method  (henceforth  abbreviated 


L-M)  is  an  approximation  b  , '  to  a  ,  of  the  form 

g(x)  g(x) 


b  =■  min  b  (n.)  , 
g(X)  Ki<k  °  1 


(4.1) 
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where 


a  "  xb  (n  )(ni  '  t0i't0i  +  ’>  ' 
a  i 


(4.2) 


with 


^Oi 


ni<1  -  "I  r  ("i  -  xi)/ni)  » 
*1-1 


(4.3) 


i.e.,  tQi  is  the  maximum  of  the  recursive  indexes  tA  defined  by  (2.4).  For  the  usual 


levels  of  a,  b  «  b^n.j).  Further,  numerical  evidence  indicates  (Harris  and  Soms 


(1983))  that  for  a  levels  of  practical  significance 

b  <  a  .  (4.4) 

g(x)  g(x) 

(4.4)  was  incorrectly  claimed  to  be  true  for  0  <  a  <  1  in  Sudakov  (1974)  and  this  is 
discussed  at  length  in  Harris  and  Soms  (1983).  However,  (4.4)  is  known  to  hold  for 
special  cases  (Winterbottom  (1974)  and  Harris  and  Soms  (1983)). 

Motivated  by  the  above,  we  now  give  an  L-M  approximation  bg(J)  to  ag(x)  for 


arbitrary  Yj^  by 


bg(i) 


min  b  (n.) 
1<i<k  °  1 


(4.5) 


where 


0  "  xb  (n  )<ni  ’  t0i't0i  +  y)  • 
a  i 


(4.6) 


with 
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j-1  3  /  j-1  3 


(4.7) 


i.e.,  t0^  is  the  maximum  of  the  recursive  indexes  t^  defined  by  (2.4).  However,  in 

this  case  it  is  not  clear  which  index  i  gives  the  minimum,  except  that  the  likely 
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candidate  is  the  one  for  which  Y ^ ,  1  <  j  <  k,  is  a  maximum.  We  might  expect,  by 


analogy,  that  for  a  levels  of  practical  interest 


bg(x)  *  ag( x ) 


5.  NUMERICAL  EXAMPLES 


For  k  =  2  and  selected  n,  Y,  x,  a  =  .05  and  .10,  Table  I  gives  ag(x) 


and  the  best  upper  bound,  u 

r  g(x) 


u  .  =  min  u  (n ,  )  , 

9<X)  1<i<k  “  1 


°  *u  (n  )*ni  “  ^0i^  *  ^0i^  +  ^ 
a  i 


and  tQi  are  defined  as  in  (4.6). 


TABLE  I. 


L-M  Approximations  and  ag(J) 


(n., ,n2 ) 

(r^Yj) 

(x1#x2) 

Qt 

bg(x) 

ag(x) 

Ug(x) 

(10,10) 

(1,2) 

(0,1) 

.05 

.3670 

.3670 

.3670 

(10,10) 

(1,2) 

(0,1) 

.10 

.4398 

.4398 

.4398 

( 10,10) 

(1,2) 

(1,1) 

.05 

.3045 

.3514 

.3670 

(  10,10) 

(1,2) 

(1,1) 

.10 

.3715 

.4227 

.4398 

(10,10) 

(1,2) 

(2,1) 

.05 

.2484 

.3151 

.3670 

(  10,10) 

(1,2) 

(2,1) 

.10 

.3088 

.3825 

.4398 

(10,15) 

(2,3) 

(0,1) 

.05 

.3695 

.3719 

.3742 

(10,15) 

(2,3) 

(0,1) 

.10 

.4425 

.4446 

.4467 

( 10,15) 

(2,3) 

(1,1) 

.05 

.2554 

.3042 

.3670 

( 10,15) 

(2,3) 

(1,1) 

.10 

.3167 

.3705 

.4398 

(10,15) 

(2,3) 

(2,1) 

.05 

.1712 

.1981 

.2431 

(10,15) 

(2,3) 

(2,0) 

.  10 

.2203 

.2513 

.3029 
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rr^ri  rsrp  >  j, j,j. 


Note  that  for  all  the  cases  in  Table  I,  b  is  a  lower  bound  for  a  The 

gi  x )  gix ) 


computations  were  done  by  a  short  FORTRAN  program,  a  listing  of  which  can  be  obtained  from 
the  author. 


CONCLUDING  REMARKS 


In  this  paper  we  have  extended  the  L-M  method  to  series  systems  with  repeated 
components.  More  work  is  needed  to  ascertain  the  region  of  validity  of  (4.8). 
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